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Abstract  A  new  test  method  for  performing  dynamic 
short-beam  shear  tests  using  a  momentum  trapped  Hop- 
kinson  pressure  bar  is  proposed.  Angle-interlock  3D 
woven  composite  specimens  were  tested  under  quasi¬ 
static  and  dynamic  loading  conditions  to  determine  the 
effect  of  loading  rate  on  damage  evolution.  The  equi¬ 
librium  condition  in  the  composite  specimen  under  dy¬ 
namic  loads  was  verified  using  finite  element  analysis 
of  the  experiment.  A  high  speed  camera  was  used  to 
capture  delamination  initiation  and  propagation  during 
both  quasi-static  and  dynamic  experiments.  Analysis  of 
the  load-deflection  curves  and  the  high  speed  images 
revealed  a  good  correlation  between  the  modes  of  dam¬ 
age  initiation  and  propagation  with  the  features  in  the 
loading  response.  The  apparent  inter-laminar  shear  strength 
and  the  bending  stiffness  increased  with  rate  of  loading. 
While  the  damage  was  observed  to  propagate  at  a 
relative  steady  rate  during  quasi-static  loading,  the  high 
rate  of  energy  input  during  dynamic  loading  resulted  in 
a  rapid  propagation  of  damage  and  a  subsequent  loss  of 
stiffness  in  the  composite  as  noted  in  the  load-deflection 
curve. 

Keywords  High  strain  rate  testing  •  Hopkinson  pressure  bar  ■ 
3D  woven  composites  •  Interlaminar  shear  strength  -  Short 
beam  shear 


T. R.  Walter  ■  G.  Subhash  (E3)  ■  B.V.  Sankar  •  M.C.  Song 
Mechanical  and  Aerospace  Engineering,  University  of  Florida, 
Gainesville,  FL  32611,  USA 

e-mail:  subhash@ufl.edu 

C.F.  Yen 

U. S.  Anny  Research  Laboratory-RDRL-WMM-B, 

Aberdeen  Proving  Ground, 

Aberdeen,  MD  21005-5069,  USA 


Introduction 

The  use  of  laminated  composite  materials  in  the  design  of 
structures  has  increased  in  recent  years  due  to  their  high 
strength  to  weight  ratio  and  increased  performance.  Howev¬ 
er,  laminated  composite  plates  often  have  low-strength  in 
the  thickness  direction  relative  to  the  in-plane  directions 
which  leads  to  a  high  propensity  for  damage  and  interlam¬ 
inar  failure.  In  particular,  impact  damage  produces  high 
inter-laminar  stresses  which  induce  delamination  damage 
[1,  2],  This  damage  may  easily  propagate  further  as  the 
composite  structure  is  exposed  to  normal  (quasi-static)  op¬ 
erating  conditions  [3,  4],  This  behavior  has  been  well  docu¬ 
mented  and  many  approaches  have  been  used  to  increase  the 
delamination  resistance  of  composites.  These  methods  in¬ 
clude  the  addition  of  through-thickness  reinforcements  such 
as  z-pinning  [5,  6],  3D  weaving  [1,  7],  toughened  epoxies 
[8-10],  and  nano-reinforcements  [11,  12],  To  determine  the 
effectiveness  of  these  techniques  on  inter-laminar  strength, 
several  test  methods  have  been  utilized.  The  double  cantile¬ 
ver  beam  (DCB)  test  [13]  has  been  used  to  investigate  the 
mode  I  fracture  toughness  of  composites  whereas  the  end- 
notch  flexure  (ENF)  test  [14]  has  been  used  to  evaluate 
mode  II  fracture  toughness.  Both  methods  require  a  pre¬ 
crack  either  produced  during  manufacture  by  inserting  a 
Teflon®  tape  or  by  cutting  a  crack  into  the  edge  of  a  sample. 
When  performing  DCB  tests  on  high  strength  composites 
(such  as  stitched  or  3D  woven  composites),  premature  fail¬ 
ure  often  occurs  due  to  the  compressive  stress  caused  by 
bending.  To  ensure  proper  failure,  a  test  method  has  been 
developed  [15]  in  which  a  combination  of  axial  and  trans¬ 
verse  loads  are  applied  to  the  test  specimen.  The  axial  tensile 
force  prohibits  the  compressive  bending  failure  while  the 
transverse  load  produces  mode  I  fracture.  Other  methods 
have  also  been  developed  to  produce  mixed-mode  loading 
[16].  The  short  beam  shear  (SBS)  test  was  developed  to 
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determine  the  inter-laminar  shear  strength  of  composites 
[17].  Unlike  the  previous  methods,  the  SBS  test  does  not 
require  a  pre-crack.  To  perform  a  short  beam  shear  test,  a 
relatively  thick  specimen  is  subjected  to  three-point  bending 
[7],  The  large  specimen  thickness  relative  to  the  test  span 
creates  a  substantial  inter-laminar  shear  stress  which  results 
in  delamination  initiation  and  propagation.  The  peak  stress 
is  calculated  assuming  a  parabolic  stress  distribution 
through  the  thickness  and  is  used  to  determine  the  inter¬ 
laminar  shear  strength  (1LSS)  of  the  composite.  However, 
the  large  contact  stresses  between  the  loading  point  and  the 
composite  often  cause  excessive  crushing  at  the  contact 
points  leading  to  premature  failure  of  the  composite,  while 
the  simplified  stress  distribution  assumption  often  misrepre¬ 
sents  the  true  ILSS.  For  this  reason,  the  strength  determined 
from  SBS  tests  are  often  referred  to  as  the  apparent  inter¬ 
laminar  shear  strength  and  can  be  used  in  the  assessment  of 
relative  perfonnance  potential  between  composites  rather 
than  as  a  material  property. 

Since  delamination  is  a  primary  mode  of  failure  in  com¬ 
posites  during  impact  or  impulse  loading,  it  is  important  to 
characterize  the  response  under  dynamic  loading  conditions. 
Several  techniques  exist  to  perform  low  velocity  impact 
tests.  The  most  common  technique  is  the  drop  tower  test 
where  a  given  impact  energy  is  imparted  to  a  composite 
specimen.  A  variety  of  composite  structures  including 
stitched  and  3D  woven  composites  have  been  tested  using 
this  method  [18-20]  at  various  impact  velocities.  Due  to  the 
combination  of  velocity  and  mass  to  impart  desired  impact 
energy,  the  rate  of  loading  is  often  non-uniform  throughout 
the  impact  event.  Also,  the  total  displacement  of  the  speci¬ 
men  is  not  controlled  and  consequently  often  varies  for  each 
test.  In  some  cases  these  qualities  may  be  undesirable  and 
therefore  alternative  methods  may  be  required.  Damage  in 
composites  can  be  evaluated  using  non-destructive  techni¬ 
ques  including  ultrasonic  measurements.  Additionally,  the 
effect  of  induced  damage  during  impact  tests  are  often 
evaluated  using  post  mortem  quasi-static  compression  after 
impact  (CAI)  tests  [21,  22].  Another  popular  method  to 
conduct  dynamic  impact  tests  on  composites  uses  gas  guns 
for  the  determination  of  impact  damage  resistance  [23,  24], 
Similar  to  drop  tower  tests,  ballistic  impact  relies  on  projec¬ 
tile  mass  and  velocity  to  impart  a  certain  level  of  impact 
energy.  In  both  of  these  tests  it  is  difficult  to  control  or 
predict  displacement  during  impact.  A  third  method  used 
to  perfonn  dynamic  tests  on  composites  is  the  split  Hopkin- 
son  pressure  bar  (SHPB)  which  is  commonly  used  to  eval¬ 
uate  high  strain  rate  response  of  many  different  engineering 
materials  [25,  26].  Traditional  SHPB  test  utilizes  elastic 
longitudinal  stress  waves  traveling  in  long  high  strength 
steel  bars  to  perform  high  strain  rate  tests.  The  sample 
shapes  may  be  varied  to  perform  a  variety  of  different  tests. 
The  SHPB  has  been  used  to  determine  the  dynamic 


compressive  strength  [27,  28],  tensile  strength  [29],  and 
shear  strength  [30,  31]  of  composites. 

Although  numerous  methods  exist  for  delamination  test¬ 
ing  of  composites,  few  methods  exist  specifically  for  deter¬ 
mination  of  delamination  strength  or  inter-laminar  shear 
strength  under  dynamic  loads.  Such  methods  are  necessary 
to  understand  the  dynamic  response  of  composites  under 
impact  loads.  In  this  manuscript  a  novel  test  method  is 
proposed  to  determine  the  apparent  delamination  strength 
of  composites  subjected  to  high  rate  of  loading  with  a  time 
scale  on  the  order  of  few  hundred  microseconds.  Experi¬ 
ments  under  static  (low  strain  rate)  loads  were  also  con¬ 
ducted  for  comparison  of  delamination  strength  and  induced 
fracture.  The  proposed  technique  demonstrates  an  alterna¬ 
tive  to  drop  tower  testing  by  producing  repeatable  deforma¬ 
tions  at  relatively  constant  velocities.  This  test  method 
allows  for  a  better  comparison  of  quasi-static  and  dynamic 
results.  A  finite  element  study  was  conducted  to  verify 
several  assumptions  used  in  the  evaluation  of  the  tests 
results. 


Experimental  Procedure 

To  evaluate  the  effect  of  loading  rate  on  the  inter-laminar 
shear  strength  of  a  3D  woven  composite,  short-beam  shear 
tests  were  performed  at  several  rates  of  loading  using  both  a 
universal  testing  machine  and  a  modified  Hopkinson  pres¬ 
sure  bar  (mHPB)  apparatus.  The  test  specimens  were  cut 
from  composite  panels  supplied  by  the  US  Army  Research 
Laboratory,  Aberdeen  Proving  Ground,  MD.  A  high  speed 
camera  was  used  during  both  quasi-static  and  dynamic  tests 
to  capture  the  deformation  and  the  modes  of  damage  in  the 
composite. 

Materials 

Angle-interlock  3D  woven  composite  panels,  manufactured 
by  TEAM  Inc.,  Woonsocket,  RI,  were  used  in  this  investi¬ 
gation.  The  preforms  were  woven  from  S-2  glass  fiber  and 
consolidated  with  SC- 15  epoxy  using  vacuum  assisted  resin 
transfer  molding.  The  weaving  pattern  consisted  of  10  fill 
tows  and  9  warp  weavers.  Each  warp  weaver  is  woven 
around  two  layers  of  fill  tows  and  the  weaving  pattern  is 
shifted  by  one  column  of  fill  for  each  row  of  weavers.  The 
pattern  repeats  every  four  rows  of  weavers  and  results  in  19 
interwoven  layers  of  S-2  glass/SCI 5  composite.  Using  rule 
of  mixtures  the  total  fiber  volume  was  estimated  to  be  48  % 
of  the  total  composite  volume.  The  architectural  details  as 
well  as  an  optical  micrograph  are  provided  in  Fig.  1(a).  The 
original  composite  panels  measured  660  mm  x  700  mm  and 
were  nominally  12.7  mm  thick.  The  as-received  panels  were 
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Fig.  1  (a)  Schematic  and  optical 
micrograph  of  3D  woven  angle 
interlock  composite  specimen, 
and  (b)  Schematic  of  static  short 
beam  shear  test  fixture.  The 
shear  stress  distribution  along 
the  center  line  is  also  shown 


sectioned  into  smaller  150  mm  x  150  mm  tiles.  These  tiles 
were  then  cut  using  a  high  speed  diamond  sectioning  saw 
into  the  desired  test  coupons  with  approximate  dimensions 
of  50  mmx20  mm.  The  20  mm  dimension  was  chosen  to 
include  at  least  two  unit  cells  in  the  width  of  the  specimen. 
The  exact  dimensions  of  each  specimen  varied  slightly  due 
to  variability  in  the  manufacturing  process.  To  remove  dam¬ 
age  induced  during  sectioning,  the  edges  of  the  specimens 
were  polished  using  standard  metallographic  techniques. 

Static  Short  Beam  Shear  Test 

The  short  beam  shear  [17]  (SBS)  test  method  was  used  to 
evaluate  the  apparent  inter-laminar  shear  strength  (ILSS)  of 
the  3D  woven  composite  at  different  rates  of  loading.  This 
test  method,  shown  in  Fig.  1(b),  uses  a  three-point  loading 
fixture  to  perform  bending  tests  on  a  test  specimen  with  a 
relatively  small  length  to  thickness  ratio.  This  geometry 
creates  a  large  inter-laminar  shear  stress  within  the  specimen 
on  either  side  of  the  central  load  point  and  promotes  delam¬ 
ination  damage.  The  inter-laminar  shear  strength  can  be 
approximated  using  the  following  formula  [17]: 

ILSS  =  0.75PB/bd  (1) 

where  PB  is  the  peak  bending  load,  b  is  the  width,  and  d  is 
the  thickness  of  the  specimen.  Quasi-static  SBS  tests  were 
first  performed  on  several  specimens  at  various  rates  of 
loading  using  a  MTS®  servo  hydraulic  universal  testing 
machine.  Figure  1(b)  shows  a  schematic  of  the  SBS  test 
fixture  with  a  fixed  support  span  of  35  mm.  Specimens  were 
tested  at  three  different  rates  of  loading;  0.0001  m/s, 
0.001  m/s,  and  0.01  m/s.  A  111  kN  load  cell  was  used  to 


measure  the  force  applied  to  each  specimen  while  the  de¬ 
flection  was  measured  using  the  LVDT  of  the  MTS®  test 
machine.  A  preload  of  200  N  was  used  to  remove  any  slack 
between  the  fixtures  and  the  specimen. 

Dynamic  Short  Beam  Shear  Test 

A  modified  Flopkinson  pressure  bar  apparatus  was  used  to 
perform  dynamic  short  beam  shear  (dSBS)  tests.  This  appa¬ 
ratus  consists  of  three  major  components;  a  single  stage  gas 
gun,  a  38. 1  mm  diameter  incident  bar  with  a  momentum  trap 
[25,  26,  32]  at  one  end  and  an  aluminum  indenter  at  the 
other  end,  and  a  rigid  anvil  which  serves  to  support  the 
specimen  and  the  force  sensors.  The  bar  material  is  a  pre¬ 
cipitation  hardened  steel  which  has  high  yield  strength,  even 
in  an  annealed  state.  A  schematic  of  the  test  system  is  shown 
in  Fig.  2(a).  The  composite  test  specimens  were  supported 
on  two  19  mm  half-cylindrical  support  rods  spanning  ap¬ 
proximately  38  mm  and  clamped  to  the  rigid  anvil.  Two 
PCB®  210b  dynamic  force  sensors  are  placed  between  the 
anvil  and  the  supports  (one  for  each  support).  These  force 
sensors  can  measure  a  combined  maximum  load  of  44.5  kN 
with  a  frequency  response  of  75  kHz.  The  sensors  are 
conditioned  using  two  Kistler  5010  charge  amplifiers  and 
are  clamped  in  place  with  a  preload  of  2.2  kN.  A  45°  wedge 
shaped  indenter  with  a  tip  radius  of  2.4  mm  is  threaded  into 
the  end  of  the  incident  bar.  The  indenter  was  made  from 
7068  Aluminum  to  reduce  the  indenter  mass,  while  main¬ 
taining  sufficient  stiffness  and  strength  to  resist  forces  dur¬ 
ing  loading.  A  350  12  foil  resistive  strain  gage  was  bonded 
near  the  midspan  of  the  incident  bar  and  was  used  to 
monitor  shape,  amplitude  and  duration  of  the  incident  wave. 


Fig.  2  (a)  Dynamic  short  beam 
shear  test  fixture  using  a  modi¬ 
fied  Hopkinson  bar  apparatus. 
Detail  of  supports,  load  cell 
placement  and  accelerometer  are 
shown  in  the  schematic  and  an 
enlarged  photo  above,  (b)  A 
schematic  is  shown  detailing  the 
process  of  the  momentum  trap 
technique 
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The  strain  gage  signal  was  conditioned  and  amplified  using 
a  Vishay  2310b  dynamic  strain  conditioner.  A  pulse  shaping 
technique  [25,  26,  33]  was  implemented  by  placing  a 
0.8  mm  thick,  19  mm  diameter  copper  disk  between  the 
striker  bar  and  the  incident  bar  before  each  test.  This  tech¬ 
nique  results  in  a  long  rise  time  in  the  incident  pulse  which 
decreases  acceleration  of  the  bar  and  reduces  signal  noise 
from  the  strain  gage  and  the  force  sensors.  Although  the 
acceleration  during  impact  is  reduced  due  to  pulse  shaping, 
the  magnitude  may  still  be  in  excess  of  200,000  m/s2.  To 
measure  this  high  acceleration,  a  PCB®  350B21  ICP  shock 
accelerometer  was  attached  to  the  end  of  the  incident  bar 
next  to  the  indenter  tip  (see  the  schematic  in  Fig.  2(a)).  The 
accelerometer  has  a  maximum  capacity  of  980,000  m/s'  and 
a  resonance  frequency  of  over  200  kHz.  This  sensor  was 
conditioned  using  a  PCB®  482A16  signal  conditioner.  Data 
was  collected  from  the  strain  gage,  force  sensors,  and  accel¬ 
erometer  using  a  LDS-Nicolet  Sigma  90-8  oscilloscope 
with  a  sample  rate  of  5  MS/s 

To  conduct  the  dSBS  test,  a  600  mm  striker  bar  was 
launched  from  the  gas  gun  toward  the  incident  bar.  Upon 
impact,  an  elastic  uniaxial  compressive  stress  wave  was 
generated  in  the  incident  bar  with  an  amplitude  and  duration 
proportional  to  the  velocity  of  impact  and  the  length  of  the 
striker  bar,  respectively.  This  compression  wave  then  travels 
down  the  length  of  the  incident  bar  until  it  reaches  the 
indenter/specimen  interface.  The  compressive  wave  causes 
the  flexure  of  the  composite  and,  due  to  the  low  impedance 
of  the  specimen,  reflects  back  as  a  tensile  wave.  In  a  tradi¬ 
tional  Hopkinson  bar,  once  the  loading  is  complete  the 
reflected  tensile  stress  pulse  travels  back  toward  the 
striker-incident  bar  interface  and  reflects  back  again  into 
the  incident  bar  as  a  compressive  wave,  which  travels  to¬ 
ward  the  indenter  and  causes  a  second  impact  on  the  spec¬ 
imen.  This  is  often  an  undesirable  consequence  of  the 
traditional  Hopkinson  bar  testing. 

If  specimens  were  to  be  recovered  after  the  test  and 
analyzed  either  through  microscopy  or  further  experimental 
testing,  it  is  important  that  the  damage  induced  during 
testing  is  limited  to  a  single  impact.  To  prevent  repeated 
loading,  a  momentum  trapping  technique  has  been  devel¬ 
oped  known  as  recovery  Hopkinson  pressure  bar  technique 
[26,  32].  This  technique  has  been  used  to  induce  controlled 
damage  in  brittle  materials  [34]  and  to  determine  dynamic 
hardness  of  several  ceramics  [25,  35,  36],  A  variation  of  this 
technique  was  used  in  this  research.  The  recovery  Hopkin¬ 
son  bar  uses  a  ‘momentum-trap’  which  consists  of  three 
components,  a  transfer  flange,  a  sleeve,  and  a  rigid  mass. 
The  transfer  flange  is  either  formed  or  machined  onto  the 
end  of  the  incident  bar.  The  flange  and  sleeve  are  impedance 
matched  to  the  incident  bar  while  the  rigid  mass  is  required 
to  be  large  enough  to  act  as  a  rigid  surface.  The  technique 
used  in  this  research  uses  a  controlled  gap  placed  between 


the  flange  and  the  sleeve.  The  details  of  this  technique  and 
the  operational  principle  are  described  below  and  are  shown 
as  a  schematic  in  Fig.  2(b). 

1)  Upon  the  striker  bar  impact  a  compression  wave  is 
generated  in  the  incident  bar  and  at  the  end  of  the  pulse 
duration,  the  gap  between  the  flange  and  the  sleeve  is 
closed  which  then  engages  the  momentum  trap. 

2)  The  pulse  reaches  the  indenter  and  after  loading  the 
specimen  reflects  back  as  a  tensile  wave  and  returns  to 
the  striker  end  of  the  incident  bar.  The  stress  wave  is 
then  transferred  by  means  of  the  flange  into  the  sleeve 
as  a  compressive  wave. 

3)  Upon  reaching  the  rigid  mass  the  stress  wave  reflects 
back  into  the  sleeve  as  a  compressive  wave  due  to  the 
large  impedance  mismatch  between  the  sleeve  and  the 
rigid  mass.  This  compressive  wave  then  returns  to  the 
transfer  flange. 

4)  Upon  reaching  the  flange  the  compressive  wave  trans¬ 
fers  into  the  incident  bar  as  a  tensile  wave  which  then 
travels  toward  the  indenter  (or  the  specimen)  and  sub¬ 
sequently  retracts  the  indenter  from  the  composite 
specimen. 

The  wave  then  repeatedly  reflects  back  and  forth  in  the 
incident  bar  such  that  it  is  always  tensile  while  traveling 
towards  the  specimen  so  that  the  incident  bar  is  incrementally 
retracted  causing  the  indenter  to  move  away  from  the  sample. 
The  result  is  a  single  dynamic  loading  on  the  composite 
specimen  followed  by  the  retraction  of  the  incident  bar. 

In  this  paper  three  possible  methods  to  determine  the 
displacement  of  the  wedge  indenter  were  examined:  (i)  a 
shock  accelerometer  mounted  on  the  incident  bar  next  to  the 
indenter,  (ii)  a  strain  gage  bonded  on  the  incident  bar  and 
(iii)  high  speed  imaging  of  the  indenter  motion  during  the 
test.  The  shock  accelerometer  measures  the  axial  accelera¬ 
tion,  from  which  velocity  and  displacement  of  the  indenter 
were  determined  by  successive  integration  of  the  accelera¬ 
tion  signal.  Integration  naturally  smoothes  the  noisy  accel¬ 
eration  signal  to  produce  smooth  velocity  and  displacement 
profiles.  The  second  method  relies  on  the  measurement  of 
the  stress  wave  using  the  strain  gage  on  the  incident  bar.  The 
strain  measured  at  a  point  on  the  bar  is  related  to  the  particle 
velocity  at  that  cross  section  through  the  following  relation: 

Up  =  Ce  (2) 

Were  e  is  the  recorded  strain  profile  of  the  incident  pulse 
and  C  is  the  longitudinal  wave  speed  in  the  bar  (-5000  m/s) 
which  may  be  determined  by  analyzing  the  time  between 
two  successive  wave  reflections  and  the  distance  the  wave 
travels  in  the  incident  bar.  The  use  of  steel  for  the  bar 
material  reduces  the  effects  of  the  dispersion  on  the  inci¬ 
dent/reflected  pulses.  In  addition  the  use  of  pulse  shaping 
reduces  the  high  frequency  components  of  the  stress  wave 
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related  to  Pochhammer-Chree  oscillations  which  further 
reduces  dispersion  and  therefore  these  effects  may  be 
neglected.  Using  this  assumption  and  the  above  relation,  it 
can  be  shown  that  the  velocity  at  the  end  of  the  incident  bar 
(i.e.,  velocity  of  the  indenter)  may  be  determined  using  the 
incident  ( Ej )  and  reflected  (er)  strain  signals  measured  at  the 
midspan  of  the  incident  bar  using  the  following  equation  [37]: 

V  =  C(£/  -  eR)  (3) 

Once  the  velocity  is  determined  the  signal  may  be  integrat¬ 
ed  to  determine  the  displacement  at  the  end  of  the  incident  bar. 
For  dSBS  tests,  the  impedance  of  the  specimen  is  negligibly 
small  compared  to  the  impedance  of  the  steel  incident  bar. 
This  results  in  very  little  loss  in  energy  during  the  loading  of 
the  specimen  and  therefore  the  reflected  stress  wave  is  almost 
equal  in  amplitude  to  the  incident  wave  but  opposite  in  sign 
(i.e.,  Er  =  —Ei ).  Using  this  assumption  we  may  simplify 
(equation  (3))  to  estimate  the  velocity  at  the  end  of  the  bar  as 

V  =  2C  Ei  (4) 

This  equation  is  useful  to  quickly  determine  the  velocity 
of  the  impact  loading  and  may  also  be  integrated  to  find  the 
displacements  at  the  end  of  the  incident  bar  for  comparison 
with  other  methods. 

Finally,  we  have  used  a  high  speed  digital  camera  (Phan¬ 
tom  v710,  Vision  Research,  Wayne,  NJ)  to  image  the  spec¬ 
imen  defonnation  during  the  impact.  The  camera  control 
software  has  the  capacity  to  manually  track  pixels  of  images 
taken  during  the  test.  Using  the  measured  thickness  of  the 
specimen  the  images  are  calibrated  to  the  proper  scale.  The 
displacement  of  the  indenter  tip  during  the  dynamic  bend 
test  can  then  be  determined  using  the  images  at  different 
time  intervals.  The  minimum  measurable  displacement  is 
limited  by  the  resolution  of  the  images  and  the  scaling  of  the 
region  of  interest.  In  this  case  the  resolved  displacement  is 
limited  to  0.108  mm  increments.  However,  the  limited  num¬ 
ber  of  data  points,  and  the  small  relative  displacements 
results  in  significant  quantization  error,  making  it  difficult 
to  accurately  determine  the  velocity  profile  from  the  high 
speed  images. 


Recall  that  the  load  applied  by  the  wedge  indenter  is 
measured  by  two  force  sensors  beneath  the  two  supports. 
The  force  measured  by  each  sensor  is  half  the  applied  load 
during  dSBS  tests.  The  signals  from  the  two  force  sensors  are 
summed  to  obtain  the  load-time  curve.  The  displacement-time 
profiles  (either  from  strain  gages,  accelerometer  or  high  speed 
images)  are  then  matched  to  the  load-time  profile  to  obtain  the 
load-displacement  curve  for  each  dSBS  test.  From  these 
curves  the  1LSS  of  the  composite  specimens  subjected  to  a 
single  dynamic  loading  is  detennined  using  (equation  (1)). 

High  Speed  Photography 

During  the  quasi-static  and  dynamic  SBS  tests,  high  speed 
imaging  was  used  to  record  damage  modes  and  damage  evo¬ 
lution  in  the  specimen.  The  camera  was  configured  with  a 
frame  size  of  1024x400  pixels  and  a  frame  rate  of  18,000  fps 
during  the  static  testing.  The  frame  size  was  reduced  for  the 
dSBS  tests  to  256x376  pixels  which  allowed  the  frame  rate  to 
increase  to  54,000  fps.  During  both  cases  a  Carl  Zeiss  Makro 
planer  100  mm  lens  was  used  with  the  aperture  set  to  f/8. 


Finite  Element  Modeling 

Equation  (1)  used  in  this  study  in  conjunction  with  the 
measured  impact  force  for  estimating  ILSS  implicitly 
assumes  that  the  specimen  is  in  quasi-static  equilibrium 
and  any  dynamic  (inertia)  effects  are  neglected.  In  order  to 
verify  this  assumption  Finite  Element  Analysis  (FEA)  of  the 
impact  test  was  performed.  The  analysis  was  simplified  by 
utilizing  the  symmetry  of  the  test  and  therefore  only  one  half 
of  the  dSBS  specimen  was  modeled  as  shown  in  Fig.  3(a).  In 
addition,  the  specimen  was  homogenized  as  an  orthotropic 
elastic  material  with  properties  given  in  Table  1  [38].  The 
use  of  fully  elastic  model  removes  any  material  rate  depen¬ 
dency  and  focuses  the  study  on  the  effect  of  inertia  and  wave 
mechanics.  The  FEA  study  included  three  rates  of  dynamic 
impact  as  well  as  quasi-static  loading  to  determine  how 
loading  rate  affects  the  outcome  of  the  analysis. 


Fig.  3  (a)  2D  FE  model  for  SBS 
test  and  (b)  Velocity  profiles  of 
indenter  for  dynamic  loading 


(a) 


(b) 
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Table  1  Material 
properties  of  the  DSBS 
specimen  for  FEA  [38] 

Density  (kg/m3) 

2300 

E,  =E2  (GPa) 

27.5 

£3  (GPa) 

11.8 

G13=G23  (GPa) 

2.14 

G12  (GPa) 

2.9 

!'13  =  V23 

0.4 

V12 

0.11 

The  commercial  FE  software  codes,  Abaqus  Standard® 
and  Abaqus  Explicit®  were  used  for  the  quasi-static  and 
dynamic  loading,  respectively.  The  specimen  was  modeled 
using  approximately  32,000  four-node  plane  strain  elements 
with  reduced  integration  (CPE4R),  and  the  specimen  thick¬ 
ness  was  set  at  20  mm.  The  indenter  and  the  support  were 
modeled  as  rigid  bodies. 

Dynamic  loading  was  simulated  using  three  different 
impact  velocities,  11,  22  and  33  m/s.  The  maximum  impact 
velocity  of  the  experiments  corresponds  to  the  lowest  rate, 
11  m/s.  The  velocity  profdes  for  the  three  impact  simula¬ 
tions  are  depicted  in  Fig.  3(b).  The  shape  of  the  profdes  was 
chosen  such  that  they  closely  resemble  the  experimentally 
measured  profde  discussed  below.  The  duration  of  each 
profile  was  modeled  to  result  in  a  final  displacement  of 
2.5  mm.  This  displacement,  which  can  be  determined  as 
the  area  under  the  v-t  diagram,  is  consistent  with  the  exper¬ 
imental  results.  Both  the  impact  force  (reaction  force  on  the 


indenter)  and  the  support  reaction  were  monitored  for  each 
time  increment.  In  the  FEA  the  displacement  of  the  rigid 
indenter  was  controlled  by  the  velocity  profile  up  to  2.5  mm 
under  the  displacement  control.  The  time  step  for  explicit 
analysis  was  calculated  automatically  by  the  FE  program 
and  it  was  in  the  order  of  10~8  s. 

The  reaction  force  curves  determined  at  the  indenter 
(impact  force)  are  shown  in  Fig.  4(a).  The  deviation  of  the 
dynamic  loading  curves  from  the  quasi-static  results  indi¬ 
cates  the  effects  of  inertia  on  the  load  histories.  From  these 
curves  it  was  observed  that  up  to  22  m/s  the  effect  of  inertia 
may  be  neglected.  At  33  m/s  this  effect  becomes  much  more 
significant  reducing  the  accuracy  of  the  results. 

Next  the  reaction  force  at  the  support  is  compared  with 
the  impact  force  at  the  indenter  to  verify  the  quasi-static 
equilibrium.  The  results  for  the  each  of  the  dynamic  loading 
simulations  are  shown  in  Fig.  4(b-d)  (note  that  the  results  for 
quasi  static  loading  were  not  analyzed  as  it  is  in  static- 
equilibrium).  If  the  specimen  is  in  equilibrium  the  reaction 
forces  at  the  indenter  and  the  total  force  at  the  supports 
should  be  equal  for  the  geometry  analyzed.  From  Fig.  4(b) 
it  can  be  noted  that  equilibrium  is  satisfied  during  impact 
with  at  a  velocity  of  11  m/s.  At  22  m/s  impact  velocity 
(Fig.  4(c)),  the  two  forces  varied  at  the  beginning  of  impact 
event  however  after  approximately  0.5  mm  of  deflection, 
equilibrium  is  reached.  It  is  observed  from  Fig.  4(d)  that 
equilibrium  is  never  fully  satisfied  when  loaded  at  a  rate 
of  33  m/s. 


Fig.  4  (a)  Load-displacement 
curves  at  the  indenter  and  (b-d) 
comparison  of  load  curves  for 
indenter  and  support  at  various 
rates  of  loading 


(a)  (b) 


Exp  Mech 


Fig.  5  (a)  Transverse  shear 
stress  distribution  in  the  speci¬ 
men  at  a  plane  midway  between 
the  indenter  and  the  support  and 
(b)  Enlarged  view  of  the  shear 
stress  distribution  near  the  center 
of  the  cross  section  (Note:  the 
theoretical  results  for  quasistatic 
and  11  m/s  are  almost  identical) 


Next  the  through  thickness  transverse  shear  stress  profile 
was  examined  for  the  quasi-static,  1 1  m/s  and  22  m/s  results. 
Since  the  results  from  33  m/s  impact  were  previously  shown 
to  be  invalid,  due  to  a  state  of  non-equilibrium,  it  was 
neglected  from  this  analysis.  These  profiles  were  determined 
at  a  point  on  a  plane  midway  between  the  support  and 
indenter  (as  shown  in  Fig.  3).  The  purpose  of  this  compar¬ 
ison  is  to  evaluate  the  effect  of  loading  rate  on  the  response 
of  the  specimen  and  compare  the  FEA  profiles  with  the 
theoretical  profile  given  by  the  equation  r  =  VQ/ It,  where 
V  is  shear  force,  Q  is  first  moment  of  area  ,  I  is  moment  of 
inertia  and  t  is  thickness,  which  is  the  basis  for  (equation 
(1)).  The  profiles  were  compared  for  a  central  deflection  of 
1.5  mm  and  are  shown  in  Fig.  5.  The  FEA  profiles  are 
shown  as  solid  lines  while  theoretical  profiles  are  shown 
as  dashed  lines.  It  is  apparent  from  Fig.  5(a)  that  there  is 
only  a  slight  deviation  between  quasi-static  and  dynamic 
shear  stress  profiles.  This  deviation  is  shown  in  an  enlarged 
view  in  Fig  5(b).  The  theoretical  results  over  predict  the 


FEA  shear  stress  by  approximately  8  %  for  quasistatic  loading 
as  well  as  for  a  loading  rate  of  11  m/s,  and  by  11  %  for  a 
loading  rate  of  22  m/s.  As  mentioned  previously  the  actual 
shear  stress  is  often  much  different  than  what  is  predicted 
using  (equation  (1))  and  therefore  is  often  used  to  compare 
different  architectures.  These  plots  reveal  that  it  is  possible  to 
compare  results  from  static  tests  with  the  dynamic  tests  at  rates 
up  to  11  m/s.  Even  at  22  m/s  the  deviation  between  the 
quasistatic  and  dynamic  results  is  only  around  3  %.  Thus, 
the  specimen  equilibrium  can  be  maintained  under  dynamic 
loads  until  an  impact  velocity  of  1 1  m/s  and  hence  the  present 
method  (equation  (1))  of  estimating  the  ILSS  of  the  material 
could  be  considered  valid. 


Results 

Typical  signals  from  the  shock  accelerometer,  strain  gage,  and 
force  sensors  during  a  dSBS  test  are  shown  in  Fig.  6(a-c).  The 


Fig.  6  (a)  Typical  acceleration 
signal  ( bold  line  is  filtered  sig¬ 
nal),  (b)  strain  profile  at  mid¬ 
point  of  the  incident  bar,  (c) 
load-time  profile,  and  (d)  com¬ 
parison  of  velocity  and  dis¬ 
placement  profiles  calculated 
from  accelerometer  ( solid  line), 
and  strain  gage  signal  ( dashed 
line),  and  displacements  from 
high  speed  images  ( circles ) 


(a) 


(b) 
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Fig.  7  Representative  load-displacement  response  of  angle  interlock 
3D  woven  composite  for  quasi-static  and  dynamic  rates  of  loading 

load-displacement  response  of  the  composite  specimen  dur¬ 
ing  the  dynamic  impact  test  was  obtained  from  these  signals.  It 
was  observed  that  the  acceleration  data  (shown  in  Fig.  6(a)) 
contains  a  large  amount  of  noise.  The  noise  may  be  reduced 
using  a  digital  filter  to  obtain  a  smooth  acceleration  curve  as 
shown  by  a  bold  line  in  Fig.  6(a).  However,  the  successive 
integration  of  the  original  acceleration  signal  reduces  the  noise 
resulting  in  smooth  velocity  and  displacement  profiles  as 
shown  in  Fig.  6(d)  without  the  need  for  digital  filters.  Using 
the  strain  gage  signal  (Fig.  6(b))  and  (equation  (4)),  the 
velocity  profile  of  the  indenter  tip  was  calculated  and 
compared  to  the  velocity  profile  from  the  accelerometer 
signal  in  Fig.  6(d).  This  profile  was  then  integrated  to 
determine  the  displacement-time  response.  Lastly,  the 
high  speed  images  captured  during  the  test  were  prop¬ 
erly  scaled  and  used  to  determine  the  displacement  of 
the  indenter  during  loading  of  the  specimen.  Comparisons  of 
the  velocity  and  displacement  profiles  determined  from  the 


accelerometer  and  strain  gage  as  well  as  displacements 
determined  from  the  high  speed  images  are  shown  in 
Fig.  6(d).  The  three  methods  indeed  show  a  good  agreement 
in  displacement  measurements. 

The  load-time  profiles  for  five  tests  shown  in  Fig.  6(c) 
demonstrate  the  repeatability  of  the  dSBS  tests.  Using  these 
profiles  and  the  displacement-time  profiles  derived  from  the 
strain  signals  the  load-displacement  curves  were  determined 
for  each  dSBS  test.  Five  tests  were  conducted  at  each  loading 
rate  and  the  curve  which  best  represents  the  average  response  of 
the  composite  was  selected  as  the  representative  load-displace¬ 
ment  curve.  Figure  7  shows  the  representative  curves  at  four 
different  loading  velocities.  The  features  in  these  curves  may  be 
related  to  the  damage  modes  observed  in  images  captured  using 
the  high  speed  camera.  These  images  are  shown  in  Figs.  8  and  9 
which  represent  typical  damage  patterns  observed  in  each  test. 
Although  slight  differences  exist  between  each  test,  several 
features  were  repeatable  and  are  discussed  below 

For  both  quasi-static  and  dynamic  tests  the  load-dis¬ 
placement  response  in  Fig.  7  shows  a  fairly  linear  region, 
followed  by  a  gradual  decrease  in  slope  of  the  loading  curve 
until  it  reaches  a  maximum  value.  It  was  observed  that  the 
slope  of  the  linear  region  is  fairly  constant  during  quasi¬ 
static  loading  at  the  selected  loading  rates,  but  increases 
significantly  under  dynamic  loading.  High  speed  images 
taken  during  this  linear  response  for  both  loading  rates  show 
small  tensile  cracks  initiated  along  the  bottom  of  the  spec¬ 
imen  as  seen  in  Figs.  8(a)  and  9(a).  These  cracks  are  due  to 
the  tensile  stress  induced  by  the  bending  but  appeared  not  to 
affect  the  overall  response  of  the  load-deflection  curve. 

After  the  linear  region  the  load  plateaus  for  a  short 
displacement  which  is  more  pronounced  in  the  dynamic 


Fig.  8  High  speed  images  of 
quasi-static  SBS  test  at  0.01  m/s. 
Magnified  images  show  detail  of 
crack  formation  including:  (a) 
tensile  cracks  near  the  bottom 
surface,  (b)  a  band  of  delamina¬ 
tion  cracks  in  the  center,  and  (c) 
delamination  propagation 
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Fig.  9  High  speed  images  of 
dSBS  test  at  10  m/s.  Magnified 
images  show  details  of  (a)  dam¬ 
age  formation,  (b)  propagation, 
and  (c)  saturation  of  damage  be¬ 
tween  indenter  and  support 


response.  The  incidence  of  this  plateau  response  coincides 
with  the  initiation  of  delamination  micro-cracks  within  the 
composite.  The  high  speed  images  of  the  0.01  m/s  test, 
shown  in  Fig.  8(a),  reveal  additional  tensile  cracks  and 
delamination  cracks  initiating  midway  between  the  indenter 
and  the  support.  As  the  material  is  further  deformed  addi¬ 
tional  delamination  cracks  initiate  between  the  warp  wea¬ 
vers  and  fill  tows  along  an  inclined  band  extending  from  the 
indenter  to  the  support  as  shown  in  Fig.  8(b).  The  dSBS 
images,  on  the  other  hand,  show  a  column  of  delamination 
cracks  which  form  near  the  center  of  the  specimen  as  shown 
in  Fig.  9(a).  While  these  cracks  were  concentrated  on  one 
side,  delamination  cracks  were  also  observed  on  both  sides 
of  the  indenter. 

As  the  displacement  increases,  the  stiffness  continues  to 
drop  with  further  increase  in  delamination  crack  length. 
After  a  displacement  of  around  1.25  mm  the  quasi-static 
load  decreases  gradually  indicating  a  further  loss  in  stiff¬ 
ness.  During  this  loading  period  the  delamination  cracks 
which  were  formed  previously,  now  propagate  over  a  large 
distance  (shown  in  Fig.  8(c))  causing  the  loss  in  load  bearing 
capacity  as  the  stiffness  of  the  specimen  is  significantly 
reduced.  During  dynamic  tests  this  drop  in  load  is  much 
more  significant,  see  Fig.  7.  Following  the  crack  initiation, 
the  propagation  of  delamination  cracks  occurs  rapidly  as 
shown  in  Fig.  9b  and  results  in  the  sudden  drop  in  load. 
The  rapid  damage  propagation  and  drop  in  load  suggested  a 
large  decrease  in  the  stiffness  of  the  specimen. 

The  above  behavior  continues  to  about  2  mm  of  central 
deflection  at  which  point  the  load  begins  to  increase  gradu¬ 
ally  until  the  test  was  tenninated  (see  Fig.  7).  In  the  dynamic 
tests,  additional  intermittent  drops  in  load  were  observed 
during  this  regime;  after  each  drop,  the  load  increases  with 
displacement.  The  increase  in  load  could  indicate  that  the 


rate  of  stiffness  loss  due  to  delamination  has  diminished, 
possibly  due  to  some  limitation  in  damage  propagation.  As 
the  delamination  cracks  approach  the  supports  or  are  hin¬ 
dered  by  the  weaving  structure  of  the  composite  the  delam¬ 
ination  growth  will  be  reduced.  Figures  8(c)  and  9(c)  show 
delamination  damage  saturation  in  the  area  between  the 
indenter  and  the  support  region. 

From  these  results  it  was  observed  that  both  stiffness  and 
peak  load  increased  as  the  rate  of  loading  was  increased  from 
the  quasi-static  to  the  dynamic  loading  regime  (Fig.  7).  The 
severity  of  the  load  drop  following  the  peak  load  was  greater 
during  dynamic  testing  compared  to  static  loading.  With  the 
high  speed  images  we  can  deduce  the  reason  for  this  large 
drop.  In  static  tests  the  damage  propagated  at  a  relatively  slow 
rate  resulting  in  a  gradual  decline  in  the  slope  of  the  load  curve 
beyond  the  peak  load.  In  the  dynamic  tests  the  rate  of  delam¬ 
ination  damage  propagation  was  much  faster  due  to  the  higher 
amount  of  energy  imparted  to  the  specimen  during  the  loading 
process.  This  energy  is  rapidly  released  during  dynamic  load¬ 
ing  through  severe  delamination  leading  to  a  loss  in  stiffness 
of  the  specimen  and  sudden  load  drop  as  seen  in  Fig.  7. 

The  maximum  load  during  each  test  was  used  to  deter¬ 
mine  the  ILSS  (equation  (1))  and  an  average  value  was 
detennined  from  five  tests  for  each  loading  rate.  The  initial 


Table  2  Apparent  ILSS  and  initial  stiffness  detennined  from  short 
beam  shear  tests 


Loading  rate  (m/s) 

Apparent  ILSS  (MPa) 

Stiffness  (kN/mm) 

0.0001 

45.40±0.75 

14.95±0.92 

0.001 

48.11±0.55 

14.40±0.40 

0.01 

50.62  ±1.23 

15.00±0.54 

10 

59.78±  1.39 

34.50±  1.17 
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stiffness  from  the  linear  response  was  also  determined  from  the 
load-displacement  curves.  Results  from  quasi-static  and  dy¬ 
namic  tests  are  summarized  in  Table  2.  It  can  be  seen  that  the 
ILSS  increases  as  the  loading  rate  was  increased  while  the 
stiffness  of  the  specimen  only  increased  under  dynamic 
loading. 

The  above  results  clearly  point  to  the  fact  that  delamina¬ 
tion  damage  is  more  severe  under  dynamic  loads  than  under 
static  loads.  Effective  strategies  are  necessary  to  limit  de¬ 
lamination  damage  in  composites  if  they  are  to  be  success¬ 
fully  implemented  in  applications  subjected  to  impact  loads. 
It  is  also  clear  that  to  understand  the  effectiveness  of  new 
composite  designs  intended  to  reduce  delamination  damage, 
controlled  tests  such  as  the  dSBS  tests  must  be  conducted  at 
high  rates  of  loading. 

Summary 

Results  from  dynamic  short  beam  shear  tests  confirmed  that 
displacement  of  the  indenter  bar  may  be  determined  accu¬ 
rately  from  stress  wave  magnitude  measured  using  the  strain 
gage  attached  to  the  incident  bar.  All  three  methods  utilized 
here  (i.e.,  accelerometer,  strain  gage,  and  high  speed  imag¬ 
ing)  for  determining  displacement  during  the  dSBS  test 
provided  consistent  results.  The  high  speed  images  provide 
a  good  method  to  validate  displacement  results  from  either 
the  strain  gage  or  the  accelerometer.  However,  the  limited 
number  of  data  points  from  the  high  speed  camera  and  the 
relatively  small  displacements  of  the  indenter  do  not  allow 
for  an  accurate  determination  of  the  velocity  profile  of  the 
indenter.  Strain  gages  provide  the  most  affordable  and  the 
easiest  method  to  implement  while  shock  accelerometers  of 
sufficient  range  are  often  limited  in  availability  and  prone  to 
frequent  damage  during  tests. 

The  interlaminar  shear  strength  determined  from  quasi¬ 
static  and  dynamic  short  beam  shear  tests  showed  a  steady 
increase  as  the  rate  of  loading  increased.  The  use  of  mo¬ 
mentum  trapping  produced  a  single  controlled  indentation 
into  the  composite  specimen.  This  technique  allowed  for  the 
recovery  of  specimens  subjected  to  a  consistent  and  meas- 
ureable  level  of  deformation.  Although  this  technique  has 
been  shown  to  produce  repeatable  results  it  is  not  intended 
to  replace  drop  tower  tests,  as  the  latter  technique  is  vital  in 
the  analysis  of  damage  resistance.  Instead,  dynamic  inden¬ 
tation  should  be  seen  as  an  additional  method  to  further 
analyze  the  rate  dependent  response  of  composites.  This 
method  for  performing  dynamic  short  beam  shear  uses  small 
specimens  to  quickly  compare  several  different  designs  of 
composites.  The  results  from  quasi-static  and  dynamic  SBS 
tests  of  3D  woven  composites  indicate  a  change  in  both  the 
load-displacement  response  and  also  in  the  initiation  and 
propagation  of  damage  characteristics. 
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